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Critical exponents for the 3D 0(n)-symmetric model with n > 3 are estimated 
on the base of six-loop renormalisation-group (RG) expansions. Simple Pade-Borel 
technique is used for resummation of RG series and Pade approximants [L/l] are 
shown to give rather good numerical results for all calculated quantities. For large 
n, the fixed point location g c and critical exponents are also determined directly 
from six-loop expansions, without addressing to resummation procedure. Analysis 
of numbers obtained shows that resummation becomes unnecessary when n exceeds 
28 provided an accuracy about 0.01 is adopted as satisfactory for g c and critical expo- 
nents. Further, results of the calculations performed are used to estimate numerical 
accuracy of the ^-expansion. The same value, n = 28, is shown to play a role of lower 
boundary of the domain where this approximation provides high-precision estimates 
for critical exponents. 
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I. INTRODUCTION 



Field-theoretical 3D 0(n)-symmetric model with self-interaction of \ip 4 type is known 
to describe critical behaviour of many basic physical systems such as Ising (n = 1) and 
Heisenberg (n = 3) ferromagnets, superffuid Bose-liquid (n = 2), polymers (n = 0), etc. In 
70-th Nickel, Meiron, and Baker, Jr. calculated all 2-point and 4-point Feynman graphs for 
this model up to six-loop order |I| paving the way for obtaining perturbative expansions of 
unprecendented length for /3-function and critical exponents. These expansions were then 
explicitly found and used, being resummed in various manners, to estimate the stable fixed 
point coordinate and numerical values of critical exponents The values obtained are 

referred today as most accurate (canonical) numbers ||. 

Explicit expressions for RG functions and numerical estimates were presented in Refs. 
]2|,|3|||] only for n = 0, 1, 2, 3. At the same time, it is desirable to have such results for n > 3. 
They are interesting from, at least, three points of view. First, there are numerous physical 
systems with many-component order parameters and these results may be relevant to their 
critical or effective critical behaviour (see, e.g. Refs. ||[7||| ) . Second, such calculations 
would enable one to clear up where resummation procedures applied to RG series become 
unnecessary, i.e. how large are the values of n for which the theory may be thought as 
possessing a small parameter. And third, high-precision numerical estimates of critical 
exponents for n ^> 1 when compared with their counterparts given by ^-expansion would 
provide an information about the numerical accuracy of this familiar approximation scheme. 

Below, six-loop perturbative expansions for /3-function and critical exponents 7] and 
7 (7 _1 ) are calculated for arbitrary n. The fixed point coordinate g c and critical expo- 
nents are estimated on the base of Pade-Borel resummation procedure and a comparison 
of these numbers with those given by unresummed RG series and ^--expansion is made. 
The outline of the paper is as follows. In Sec. |TJ the renormalisation scheme is formulated, 
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RG expansions are written down, the resummation technique is described, and numerical 



results obtained are collected. In Sec. Ill they are discussed along with their analogues 



resulting from unresummed six-loop series and --expansion and corresponding inferences 



are presented. Section IV contains conclusions 



II. RG SERIES AND NUMERICAL RESULTS 



The Hamiltonian of the model to be studied reads: 



// - - / d'.r 



where ip a is a vector order parameter field, a = 1, . . . ,n, a bare mass squared tuq being 
proportional to the deviation from the mean-field transition point. 

We calculate the /3-function and critical exponents within a massive theory. The renor- 
malized Green function Gn(p,m, g) and four-point vertex function Tn(p,m, g) are normal- 
ized at zero momenta in a conventional way: 



9G^-(p,m,g) 
dp 2 

T R (0,m,g) = mg 
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= 1 , (2.2) 

p 2 =0 



with one extra condition being imposed on the ip 2 insertion: 



rf(p,q,m,g) p=g=Q = l ■ (2.3) 

Since combinatorial factors and momentum integrals for 2-point and 4-point Feynman 
graphs are known [TJ the calculation of the /3-function and critical exponents (anomalous 
dimensions) within six-loop approximation is straightforward (see, e.g. 0). The results are 
as follows: 



(3(g) = g- g 2 + — ^6.0740740877 + 28.14814815V - - — ^—^ ( 1.34894276n 2 
[n + 8) 2 \ ) (n + 8)^ V 

+54.94037698n + 199.6404170 V + - — ^— f-0.15564589n 3 + 35.82020378n 2 

/ (n + 8) 4 V 

+602.521230577 + 1832.206732V - - — 1 —— fo.05123618n 4 + 3.23787620n 3 

r (n + 8) 5 V 

+668.5543368n 2 + 7819.564764n + 20770.17697V + — - f-0.02342417n 5 



r (n + 8) 6 , 

-1.07179839n 4 + 265.8357032n 3 + 12669.22119n 2 + 114181.4357n 
+271300.0372^ 7 , 

rj(g) = - — 1 —— fo.2962962964n + 0.5925925928V 2 + - — — fo.0246840014n 2 
(n + 8) 2 \ ) (n + 8) 6 \ 

+0.246840014n + 0.3949440224 V + - — ^— f-0.0042985626n 3 



r (n + 8) 

+0.6679859202n 2 + 4.609221057n + 6.512109933 V 4 - - — ^— ( 0.0065509222n 4 



/ (n + 8) 1 

-0.1324510614n 3 + 1.891139282n 2 + 15.18809340n + 21.64720643V 



+ - — ^—/-0.0055489202n 5 - 0.0203994485n 4 + 3.054030987n 3 
(n + 8) 6 V 

+64.07744656n 2 + 300.7208933n + 369.7130739^ 6 , 

77+2 77 + 2 1 / 

= 1 " 0( + T^W* - -— — 0.879558892677 2 + 6.485476868n 

2(77 + 8) (77 + 8) Z (77 + 8) 6 V 

+9.452718166 V + - — 1 — 7 f-0.1283321043n 3 + 7.96674070377 2 
r (t7 + 8) 4 V 

+51.8442129877 + 70.79480631^1 g 4 - fo.0490966058n 4 

(77 + 8) 5 V 

+4.288152493n 3 + 108.3618219n 2 + 537.8136105n + 675.6996077V 



f- — 1 — 7 - ( -0.025926794577 5 - 1.61862784377 4 + 85.54569746n 3 



[n + 8) 6 

+ 1538.81823577 2 + 6653.956526n + 7862.074086^ 6 . 

These series are known to be divergent (asymptotic). To extract the physical information 
which they contain some resummation procedure should be employed. We use the Pade- 
Borel method, i.e. construct Pade approximants [L/M] for Borel transforms which are 
related to functions to be found ( "sum of series" ) by the formula 



oo 

f(x) = J2c k x k = / e-*F(xt)dt , (2.7) 



k=0 



nv) = > (2-8) 

fc=0 



and then evaluate the integral (|2.7| ) where series (|2.S| ) possessing non-zero radii of conver- 
gence are replaced by corresponding Pade approximants. 

Starting from six-loop expansions available it is possible to construct different sets of 
Pade approximants: [L/l], [L — 1/2], etc., where L = 6 for /3-function and L = 5 for critical 
exponents. As we found, approximants 

[L/l] = (l + b 1 y)- 1 j2a i y i (2.9) 

i=0 

which generate following expressions for sums of the series 

f(x) = ze- z Ei(z) ]T aii-h)-* ~ E ^{-b^ ^ k\z~ k , 

i=0 i=l fc=0 



Z = — 



z 

— , = / e'r^t (2.10) 



give the best results. They are presented in Table |. The estimates for 7 and 77 originate 
from series ( |2.5| ) and ( p.6| ) while numerical values of critical exponents u, a, and f3 were 
determined by means of well-known scaling relations. The exponent 7 was calculated also 
via resummed RG expansion for the exponent r] 2 = (1 — 7) (2 — rj)/^ and numbers were 
obtained which differ from those resulting from ( |2.6| ) by no more than 0.003; corresponding 
averages stand in Table |. This table contains as well, for comparison, numerical results 
found earlier for n = 0, 1,2,3 on the base of higher-order RG expansions in 3 and 4 — e 
dimensions using alternative resummation techniques |2|,0,|TIJ. It is worthy to discuss these 
results along with ours in more detail. 

As we can see, there are small differences between our estimates and their counterparts 
obtained in Refs. from 3D RG expansions of the same length. They are caused by 



use of different resummation procedures. Indeed, the authors of Refs. |2||| employed the 
Borel-Leroy transformation 

oo 

f(x) = J t B e~ i J : {xt)dt (2.11) 
o 

instead of Eq. ( |2.7| ) in their calculations. The parameter B was chosen to meet the known 



large-order behaviour of coefficients Ck in perturbative expansions [11,12 



c k ~ k\{-afk b , k -> oo , (2.12) 

where a = 0.147774 for the model ( ^.1| ) and b is equal to 2 + |, or 3 + ~, or 5 + | depending 
on the RG function expanded. We use much simpler method which ignores some part of 
information ( [2.12j ) but leads, nevertheless, to numerical results rather close to those given 
by more sophisticated techniques. It is not surprising since the main property of c& - their 
factorial growth, is taken into account in our analysis while the rest of information about 
Ck being incorporated enables one to reduce the apparent errors of estimation keeping the 
location of fixed point and critical exponents practically unchanged (see, e.g. Ref. || for 
detail). Dealing with simple Pade approximants [X/l], we avoid also, to a certain extent, the 
problem of poles. The point is that these approximants turn out to have no real and positive 
poles for n < 38 in the case of critical exponents and up to n = 80 for the /^-function. That 
is why Table | ends at n = 32. Since for n = 0,1,2,3 our procedure gives critical exponents 
values which are almost identical to known high-precision estimates |2|,|3|,[1(J , we believe that 
the rest of the results listed in this table are also very close to exact numbers. 

III. LARGE N AND ^ EXPANSION 

How can we estimate g c and critical exponents for n > 30 ? It is well known (and clearly 
seen from Eqs. (f2.4|)— (P75|)) that coefficients of RG expansions are decreasing when n grows 
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up. Hence, for large enough n the theory should possess a true small parameter as, say, the 
quantum electrodynamics does. In such a case, all quantities of interest can be obtained 
directly from corresponding perturbative expansions, without addressing to resummation 
technique. To find the minimal value of n which may be referred to as "large enough" we 
have calculated g c for 20 < n < 60 using original and Pade-Borel-resummed series (|2.4j ). (It 
should be reminded that the approximant [6/1] for the Borel transform of /3— function has 
no dangerous poles within this segment.) The results are presented in Table [TT]. Values of g c 
given by these two approximations are seen to differ from each other by 0.9% for n = 28 and 
this difference diminishes rapidly with increasing n. So, if the accuracy of order of 1% for g c 
was adopted as satisfactory, the resummation of six-loop expansion for /3-function becomes 
unnecessary when n exceeds 28. The same turns out to be truth for the critical exponent 7 
as is seen from Table |TTT| (the first and the second lines). 

For large n, another approximate method may be used to calculate critical exponents. 
We mean famous --expansion. Within the second order in - exponents 7 and r\ are known 



Series for other critical exponents are easily obtained via scaling relations. 

It is interesting to evaluate the accuracy of numerical results given by ^--expansion. 
We can get such an information comparing numbers resulting from Eqs. (|3.1| ) and ( |3.2| ) 
for various n with their counterparts obtained on the base of resummed (n < 32) and 
unresummed (n > 32) six-loop RG series. On the other hand, this comparison would help 
us to determine the accuracy of the employed approximation itself in the limit n —> 00 where 
--expansion's results are exact. 



to be H3: 




(3.1) 



8 1 512 1 



(3.2) 



V 



3?r 2 n 27tt 4 n 2 
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Corresponding estimates for exponent 7 are listed in Table [II. These numbers show 
that numerical accuracy of Eq. ( |3.1|) becomes better than 1% when n exceeds 28. Values of 
77 given by six-loop RG series and Eq. ( ^.2| ) are very small and not presented here. They 
differ from each other by approximately 10% for n > 28. Moreover, this discrepancy persists 
up to largest values of n studied. It is not surprise. The point is that, for extremely large 
n, only leading terms in n contribute to 77 in each order in g. Since g c = 1 + 0(n _1 ), g c 
should be put equal to unity within this limit. Hence, corresponding total contribution in 
the case of six-loop RG series may be found by summing of coefficients of all leading terms 
in Eq. (|2.5|) . Such a procedure gives 77 = 0,3 ° 458 , while the exact asymptotic expression 
resulting from Eq. ( ^.2[ ) is 77 = °- 2 ™ 19 . So, the approximate asymptotic estimate for 77 differs 
from the exact one by 13%. This difference, however, practically doesn't influence upon 
numerical values of other critical exponents calculable by scaling relations since for n > 28 
the exponent 77 < 0.01. 

We see that simple formulas ( |3.1| ) and ( |3.2j ) enable one to estimate all critical exponents 
for the model (|2.1|) with an accuracy of order of 0.01 provided n > 28. Moreover, for such 
n second-order terms in these formulas may be, in fact, neglected since their contributions 
are very small. 



IV. CONCLUSION 

Critical exponents of the 3D 0(n)-symmetric model have been estimated from six-loop 
RG series for n > 3. RG expansions have been resummed by means of simple Pade-Borel 
technique and approximants [6/1] (/3-function) and [5/1] (critical exponents) have been 
shown to provide rather good numerical results for all calculated quantities. It has been 
found that for n > 28 the theory may be thought as possessing a small parameter, i.e. the 
fixed point coordinate and critical exponents may be determined with errors about 0.01 or 
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less directly from higher-order RG series, without use of resummation procedure. Numerical 
accuracy of the ^-expansion has been also estimated. The same value, n = 28, has been 
shown to play a role of a lower boundary of the region where this approximation provides 
high-precision results for critical exponents. 
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TABLES 

TABLE I. The stable fixed point location and critical exponents obtained within six-loop 
approximation using the Pade-Borel resummation technique. 



n 


9c 


7 


V 


V 


a 


P 





1.402 


1.160 


0.034 


0.589 


0.231 


0.305 




1 491* 


1 161 * 

1.1U1 


096* 


588* 


n 9^6* 

L/. Z*OU 


S09* 




1 41 7t 


1 1fi9t 


n 09fii 

U.UZU 


^RR^ 
u.ooo 




n inot 


1 


1.401 


1.239 


0.038 


0.631 


0.107 


0.327 




1 A 1 fit 


1 9/1 1 1 


n n^i t 


U.DOU 1 


n 1 1 nt 

U. 1 1U 1 






1.414* 


1.240* 


0S2* 


6S0* 




S25* 








U.UoO 


U.DZo 






2 


1.394 


1.315 


0.039 


0.670 


- 0.010 


0.348 






1 fit 


U.UoZ 1 


u.ooy 1 


n nn7t 

- U.UU ( 1 


n ^/ifit 




1.405* 


1.316* 


0.034* 


0.669* 




0.346* 








0.037* 


0.665* 






3 


1.383 


1.386 


0.038 


0.706 


- 0.117 


0.366 




1.392* 


1.390* 


0.031* 


0.705* 


- 0.115* 


0.362* 




1.391* 


1.387* 


0.034* 


0.705* 




0.365* 



12 



0.037* 0.698* 



4 




1.369 


1.449 


0.036 


0.738 


- 0.213 




0.382 


5 




1.353 


1.506 


0.034 


0.766 


- 0.297 




0.396 


6 




1.336 


1.556 


0.031 


0.790 


- 0.370 




0.407 


7 




1.319 


1.599 


0.029 


0.811 


- 0.434 




0.417 


8 




1.303 


1.637 


0.027 


0.830 


- 0.489 




0.426 


9 




1.288 


1.669 


0.025 


0.845 


- 0.536 




0.433 


10 




1.274 


1.697 


0.024 


0.859 


- 0.576 




0.440 


12 




1.248 


1.743 


0.021 


0.881 


- 0.643 




0.450 


14 




1.226 


1.779 


0.019 


0.898 


- 0.693 




0.457 


16 




1.207 


1.807 


0.017 


0.911 


- 0.732 




0.463 


18 




1.191 


1.829 


0.015 


0.921 


- 0.764 




0.468 


20 




1.177 


1.847 


0.014 


0.930 


- 0.789 




0.471 


24 




1.154 


1.874 


0.012 


0.942 


- 0.827 




0.477 


28 




1.136 


1.893 


0.010 


0.951 


- 0.854 




0.481 


32 




1.122 


1.908 


0.009 


0.958 


- 0.875 




0.483 


t Quoted from Ref. @. 

* Quoted from Ref. @. 

* Quoted from Ref. |10| 














TABLE II. Coordinates of the fixed point obtained from Eq. 
resummation procedure (PB) and by direct summation (DS). 


pi]) with use 


of Pade 


-Borel 


n 




20 


24 


28 32 


36 


40 


50 


60 


9c 


(DS) 
(PB) 


1.2184 
1.1768 


1.1725 
1.1538 


1.1458 1.1273 
1.1359 1.1216 


1.1134 
1.1099 


1.1025 1.0830 
1.1003 1.0822 


1.0699 
1.0696 
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TABLE III. Values of the critical exponent 7 obtained by direct summation of the RG expan- 
sion (DS), by means of Pade-Borel resummation technique (PB) and from ^-expansion (Eq. ( |3.lD ). 



n 


20 


24 


28 


32 


36 


40 


50 


70 


100 


500 


(DS) 


1.8990 


1.8991 


1.9075 


1.9165 


1.9245 


1.9314 


1.9447 


1.9606 


1.9725 


1.9946 


(PB) 


1.8466 


1.8737 


1.8932 


1.9078 


1.9222 












j_ 

n 


1.8702 


1.8930 


1.9090 


1.9208 


1.9299 


1.9372 


1.9501 


1.9646 


1.9754 


1.9951 
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